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ABSTRACT 

We present an analysis of the Globular Cluster (GC) population of the elliptical 
galaxy NGC 4261 based on HST WFPC2 data in the B, V and / bands. We study 
the spatial distribution of the GCs in order to probe the anisotropy in the azimuthal 
distribution of the discrete X-ray sources in the galaxy revealed by Chandra images 
(^Zczas c t al. 2003). The luminosity function of our GC sample (complete at the 90% 
level for m^, = 23.8 mag) peaks at m^, = 25.llg g mag, which corresponds to a distance 
consistent with previous measurements. The colour distribution can be interpreted as 
being the superposition of a blue and red GC component with average colours V — I 



1.01 



-0.06 



mag and 1.27_Q Qg mag, respectively. This is consistent with a bimodal 
colour distribution typical of elliptical galaxies. The red GCs radial profile is steeper 
than that of the galaxy surface brightness, while the profile of the blue subpopulation 
looks more consistent with it. The most striking finding is the significant asymmetry in 
the azimuthal distribution of the GC population about a NE-SW direction. The lack 
of any obvious feature in the morphology of the galaxy suggests that the asymmetry 
could be the result of an interaction or a merger. 
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1 INTRODUCTION 

Elliptical galaxies are spheroidal systems which are expected 
to present a highly uniform spatial distribution of starlight 
and Globular Cluster (GC) sources. Nevertheless, optical 
images show that several elliptical galaxies exhibit struc- 
tures such as shells, ripples, arcs, tidal tails and other asym- 
metric features fe.g. ISchweizer fc Seitzerlll992f ). One of the 
most popular scenarios for elliptical galaxy formation sug- 
gests that these features are the result of major mergers. 
Num erical simulations (e.g. lToomre &: Toomreiri972l : lBarnesl 
llQSS 'l show indeed that such mergers are characterized by 
peculiar remnant structures of this type. Minor merging 
events and interactions can subsequently disturb further the 
galaxy morphology. 



Merger events and galaxy encounters affect the spa- 
tial distribution of the stellar population as well as that 
of the GC population. ISchweizer fc Seitzed (|l992l ) discov- 
ered significant structure in the starlight of several elliptical 
galaxies, and disturbed orbits are found in NGC1404 and 
NGC1399, probably d ue to a dynamical interaction between 
the two galajcies ([Napolitano. Arnaboldi. fc Capacciolil 
l2002l : iBekki et al]|2003l '). The amplitude of the deviation 
from the isotropy in either the stellar or GC populations 
in elliptical systems obviously depends on the time lapsed 
since the merging/interaction event, as well as on the na- 
ture of the event itself, and so may be used to explore the 
elliptical galaxy's formation history. 
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1.1 The Case of NGC4261 



NGC 4261 is a nearby (29.4±2.6 Mpc; IJensen et al 
200S ) early-type galaxy in the Virgo W cloud l|Garcia 



19931 : iNoltheniujfigga V It is a well-studied object thanks 
to many interesting features: it has an active nucleus 
hosting a supermassive bla ck hole (4.9±1.0) x 10*Mq 
iFerrarese. Ford, fc Jaff j Il996h . it hosts the radio source 
3C 270, whose prominent radio jets (|Birkinshaw fc DaviesI 
1 19851 ) are associated wit h X-ray emission close to the nu- 
cleus jWorrall et al.ll2010|'). It shows a 20"dust lane along the 
north-south axis (|Martel et aIll2000D and has an impressive 
nucle ar dust disk, discovered with HST WFPC2 (|jaffe et al.l 
1 19961 ). 

Apart from the "boxy" isophotes fe.g. iNieto fc Bended 

1 19891 ) and weak evidence for a tidal arm in the NW di- 
rection (|Tal et al.1 [2OO9I ). this galaxy does not show any 
definite story of recent interaction such as shells, ripples, 
rings, resulting in a low fi ne structure parameter (E=1.0; 
ISchweizer fc Seitzeij Il992i ). Therefore, despit e being the 
most massive object in a poor group of galaxies (|Davis et al.l 
1 19951 ). we can claim that NGC 4261 does not show evidence 
of recent gravitational interaction. 

Given the undisturbed starlight distribution of 
NGC 4261, the discrete X-ray sources (i.e. X-ray bina- 
ries) would be expected to be distributed uniformly. How- 
ever, Chandra data pointed to an interesting anisotropy in 
their azimuthal distribution, and more precisely to an ex - 
cess between P A. =14 0° and P.A.=190° i|Zezas et al.ll2003l ). 
I Giordano et al.l (|2005l ) associated 50% of these sources with 
GCs identified in archival optical images and suggested that 
the Chandra sources are accreting Low Mass X-ray Binaries 
(LMXBs). This result suggests a non-uniform distribution 
of the GC population. 

In order to investigate this hypothesis, we requested 
deep HST WFPC2 data (proposal ID 11339; PI: Zezas, A 
- Table [1]) to characterize the GC population of NGC 4261. 
In this paper we present a study of the radial and azimuthal 
distributions of the blue (metal poor) and red (metal rich) 
GC subpopulations. We assess the extent of the asymmetry 
and discuss its origin in the context of a galaxy interaction. 
In a following paper (Bonfini et al., in preparation), we will 
extend the analysis to the LMXB population using the new 
deep Chandra pointing. In particular, by probing whether 
the LMXBs and GCs indeed have similar spatial distribu- 
tion, we will add to the deba te on whether GCs are the 
sole b irthplaces of LMXBs (e.g. IWhite. Sarazin. fc Kulkarnil 
20021). or field LMXBs form in dependently in situ (e.g. 



Kundu. Maccarone. fc Zepj I2OO7I . and references therein). 



or, as suggested by the recent work of iKim et al.l (|2009l ). 
both formation processes contribute to the field LMXB pop- 
ulation of elliptical galaxies. 

The outline of the paper is as follows. In 33 we de- 
scribe the HST data reduction and source detection. In fj3] 
we report the criteria used to define GC candidates and we 
describe the procedure used to correct the Luminosity Func- 
tion (LF) for incompleteness. In ij4] we study the Globular 
Cluster Luminosity Function (GCLF) and colour distribu- 
tion. In 35] we study the azimuthal and radial distributions 
of the GCs. In Sj6]we discuss our results in the context of the 
recent interaction history of NGC 4261. 
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Figure 1. / band mosaic. Blue and red markers represent the 
position of blue and red secure GCs down to the 75% complete- 
ness Umit (sample and choice of completeness limit are defined in 
Section ^and JS] respectively). These are the sources on which 
we performed the study of the spatial distribution. The white 
open circles represent the lo cations of X-ray sources in the source 
list of lGiordano et ahllioosl . 



2 OPTICAL OBSERVATIONS AND DATA 
ANALYSIS 

2.1 Data Reduction 

HST WFPC2 data for NGC 4261 were obtained with the 
F450W, F606W and F814W filters (roughly corresponding 
to the B, V and I bands, respectively; see WFPC2 Instru- 
ment Handbook, Version 10.0, Table S.lQ ). The filters were 
chosen to provide the highest sensitivity while allowing us 
to distinguish the blue and red GC subpopulations. 

In order to reject cosmic rays the observations were 
performed following a two point dither pattern with a 0.3 
pixel offset. A total of 12 exposures of 400 s each were 
acquired in each filter in a 3x2 grid. The pointings were 
planned in order to cover the D25 (major diameter; 4'. 07 
Ide Vaucouleurs et al.|[l99ll ) area of the galaxy. The details 
of the observations are reported in Table [1] 

Mosaics of the HST images were created with the fol- 
lowing procedure. The task MultiDrizzle (within the PyRAF 
package, DITHER version 2.3) had been first used to pro- 
duce the distortion corrected frames for each exposure. 
Then, the relative offsets between the fields were determined 



^ |http://www.stsci.edu/hst/wfpc2/documents/IHB_17.html| 
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Table 1. Log Of The Observations 



Instrument 


Filter 


Number 
of 
Fields 


Date 


Exposure 
per Field 
(sec) 


HST WFPC2 

HST WFPC2 
HST WFPC2 


F450W (B) 

F606W (V) 
F814W {/) 


6 

6 
6 


2008 Jan 12 
2008 Jan 13 

2007 Dec 26 

2008 Jan 28 
2008 Mar 10 


{2x)400 

(2x)400 
{2x)400 



by comparing the coordinates of stars present in overlap- 
ping regions of these frames. The pointlike sources within 
each frame were detected using daofind, while the PyRAF 
tasks xyxymatch and geomap were used to calculate the off- 
sets. Due to the small number of stars in the field of view, 
we could identify only 2-5 sources in common between each 
pair of fields, therefore in some cases we had to include some 
GCs as reference objects. The offset table created in this way 
was fed to MuHiDrizzle in order to create the final version of 
the mosaic. The mosaics for the different filters were regis- 
tered with respect to the / band image. The final image was 
binned to the Wide Field Camera CCD pixel scale (WFC; 
nominally 0.0996" /pixel). 

Subsequently, we corrected for systematic astrometry 
errors by cross-correlating the HST coordinates of the 50 
brightest objects (counts > 4000) in the / band mosaic 
(see i|2.2l for the source detection technique) against the 
coordinates of the 50 brightest Sloan Digital Sky Survey 
- Da ta Release 5 (SDSS DR5; lAdelman-McCarthv et alj 
I2OO7I ) sources in the field. The SDSS catalogue pro- 
vided a more reliable calibration since it is richer in as- 
trometric reference stars than the HST field. We per- 
formed the astrometric correction using the WCSTOOLS 
imwcs, which produced a significant number of matches 
(11/50), allowing us to measure a negligible WCS offset of 
(ARA,ADec) = (-0.02±0.10",0.01±0.13"). The same WCS 
solution was then applied to the B and V band mosaics. 

Figure[T]shows the final drizzled image for filter /, along 
with the location of the blue and red GCs of our secure sam- 
ple down to the 75% completeness limit (the sample and 
choice of completeness limit are defined in 32] and ^ re- 
spectively). In the figure we also show the locati ons of the 
X-ray sources in the list of [Giordano et al] (|2005l ). 

2.2 Source Identification and Photometry 

We used SEXTRACTOR (|Bertin fc Arnoutd figg^ l to de- 
tect discrete sources in the mosaics. This package performs 
source detection and provides net counts for each source 
after e stimating and subtr acting the local background. Fol- 
lowing (^rbeseTal] l|2004l 'l. who used SEXTRACTOR for the 
detection of GCs in WFPC2 data of nearby galaxies, we chose 
the following software parameters: minimum detection area 
of 4 pixels; detection threshold of 1.5 a; threshold for pho- 
tometry analysis of 1.5 a; background grid of 16x16 pixels; 
filtering Gaussian of 2.5 pixels Full Width at Half Maximum 
(FWHM; roughly equal to the value of the FWHM of the 
point-like sources in the mosaics). 

A typical GC ha s a half fight radius rh—3 pc (e.g. 
I Ashman fc Zepj Il998l '). At the distance of NGC 4261 this 



value translates into a size of 0.02"or 0.2 pixels at the scale 
of our mosaics, meaning that the GCs in our images (in- 
cluding the wings in their light profile) are only marginally 
resolved, which does not allow accurate comparison to the 
Point Spread Function (PSF; PSF FWHM ~ 2 pixels in 
our mosaics) to determine their extent. Considering the low 
luminosity of the sources, we decided to perform fixed aper- 
ture photometry in order to minimize the photometric error. 
This decision is justified by the results of the incompleteness 
simulation (see i]3.2p . We adopted an aperture diameter of 
5 pixels (or 0. 5"), encompassing ^95 % of the total flux of a 
point source (|Holtzman et al.lll995al ). 

The conversion from instrumental magnitudes to the 
UBVRI system was perf ormed using the WFPC2 calibra- 
tions of iHoltzman et al. I ([l995a). When available, trans- 
formation parameters based on observational data were 
preferred over synthetic ones. The B-V and V-I colours 
were obtained f r om an analytic solution of equation (8) of 
iHoltzman et al.l (|l995al ) for the single B, V and / bands: 

B - V = (1) 

^'•^i.FS - ^2.FS> 



fv 

where Ti^fs and T2,fs are the Hi^ht-to-U BV RI transfor- 
mation coefficients, Zfs the zeropoints from lHoltzman et al.l 
(|l995ah and /s , fv and // correspond to the aperture fiuxes 
in the B, V and I band respectively. 

Aperture corrections (evaluated with the method de- 
scribed in H3.2p were applied to account for aperture losses. 



3 DEFINITION OF THE GLOBULAR 
CLUSTER POPULATION 

3.1 Globular Cluster Selection 

We created our GCs list by cross-correlating the 
SEXTRACTOR output catalogues for the V and I band. The 
B band catalogue was excluded from the cross-correlation 
due to the limited number of sources detected in this band 
(because of the relatively red colour of GCs). The result of 
the match was subsequently screened according to the fol- 
lowing criteria: 

Selection on position. We restricted the sample to 
sources within the D25 ellipse of the galaxy (4'. 07 major 
diameter; 0.8 axis ratio b/a). The SEXTRACTOR source de- 
tection efficiency dropped drastically at the galaxy center 
due to the high background, with virtually no detections 
within a galactocentric radius of Roc ~ 25". 

Selection on axial ratio. GCs are usually close to perfect 
spherical systems, therefore we could co nfidently exclude ob - 
jects with significant axial ratios. As in I Jordan et al.l (|2004l ) 
we adopted a mean axial ratio (g) (i.e., the average of the 
axial ratios in different filters) and we set conservative limits 
0.5 < (G) < 2.0. 
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Selection on (S/N). We limited the signal to noise ratio 
of each source to be higher than 20 for both the V and / 
bantfl. 

Selection on FWHM. As mentioned in Section 12.21 
GCs at the distance of NGC 4261 are expected to be only 
marginally resolved. Therefore, the minimum FWHM is dic- 
tated by the PSF FWHM of the EST WFPC2 cam era, which 
is around 2-2.5 pixels fsee iHoltzman et al.lll995bl . Figure 5). 
Quantifying the upper limit, instead, is a more difficult task 
since it should take into account the convolution of the ra- 
dial profile of an extended GC with the instrument PSF. 
Based on extensive simulations (see i]3.2p . we set the upper 
limit to a few times the PSF FWHM in order to account 
for the largest GCs. Larger objects are expected to be back- 
ground galaxies and therefore must be rejected. Based on 
these considerations, we restricted the FWHM within the 
limits: 1.5<FWHM<4.0 pixels, for every band. 

Selection on colour. The colour was constrained to be 
in th e ran ge observed in elliptical gal axy GCs: 0.6<y- 7<1.6 
(e.g. [Ashm an & Zcpi 1998; Kundu fc WhitmordlioOll '). 

The forementioned criteria resulted in a sample of 
718 "secure GCs" (down to a minimum V magnitude of 
~25.4 mag). A second, less strict sample, was defined from 
the SEXTRACTOR catalogue for the V band applying the 
same selection except for the colour criterion. In this way, we 
found a total of 1067 "GC candidates" (down to a minimum 
V magnitude of ~25.8 mag). 

False detections associated with background galaxies, 
although not excluded, are expected to be minimal, due to 
the restrictions applied on the axis ratio and FWHM and, 
in the case of the secure sample, to the colour limits (late 
type galaxies are much bluer than the typical GC). As a 
further test, we investigated the effectiveness of the FWHM 
selection in rejecting background galaxies as opposed to an 
equivalent criterium often used in similar HST studies of ex- 
tragalactic GCs, which is the difference in magnitude within 
two different apertures. To do so, we performed a compar- 
ison of the flux ratio between a 1.5 pixels and a 2.5 pix- 
els aperture (roughly the size of the PSF FWHM) against 
the FWHM of the candidate GCs sample. This comparison 
showed a tight correlation between the two quantities (al- 
though affected by significant scatter). Therefore, applying 
either of the two criteria would lead to the same results. We 
also verified that source confusion (i.e. blending of sources) 
is negligible, if present at all. 

Table |4] lists the photometric parameters of the objects 
in the secure GCs sample. When available, the B magnitude 
and the B-V colour have been included. 



3.2 Evaluating Incompleteness 

In order to evaluate the incompleteness of the GC samples 
(i.e., the fraction of GCs not detected due to faintness or is- 
sues related to the detection process) , we set up an artificial 
source test to calibrate the SEXTRACTOR results. Simu- 
lated GCs were added to the NGC 4261 HST mosaics and 



2 The SEXTRACTOR thresholds discussed above refer to the 
(S/N) for any individual pixel of a source, rather than to the 
(S/N) for the total flux within the aperture. 



their characteristics were measured with SEXTRACTOR us- 
ing the same setup as for the real data. The simulation was 
repeated several times in order to improve the statistical re- 
sults. The details of the artificial source test are described 
in Appendix 1X1 

Figure [2] reports the simulation results for the V (top) 
and / (bottom) bands. The left panels of the figure show 
the detection rate (i.e. percentage of detected sources over 
simulated sources) as a function of the input magnitude. 
Input sources are detected with an efficiency of ~90% (or 
higher) down to 23.8 mag and 22.7 mag for the V and / 
band respectively. We define these magnitudes as the de- 
tection thresholds. In order to assess the dependence of the 
completeness curves on the background light we repeated 
the simulation within 3 annulii located at different galac- 
tocentric radii (between the simulation limits Rgc ~ 25" 
and Rgc ~ ^25). We found that the detection efficiency 
indeed shows a dependence on the brightness of the back- 
ground, although the differences are significant only well be- 
low the 80% completeness level. The completeness correction 
applied to the V band histogram in ii4.1l took this effect into 
account. Since the detection efficiency drops drastically for 
Rgc < 25 ", we adopted this as the minimum galactocentric 
radius for the study of the spatial distribution of the sources 
(3Sll- The right panels of Figure[2]show the ratio of the 0.5" 
diameter aperture fiux measured with SEXTRACTOR over 
the total input fiux of the simulated sources, as a function 
of the input magnitude. At the faintest magnitudes, the fiux 
ratio deviates from the constant value since the Poissonian 
variation of the background counts becomes an important 
contribution to the signal within the aperture. For this rea- 
son, we estimated the aperture correction only down to the 
75% completeness limit. We find that our fixed 0.5" aper- 
ture encompasses 79%, 78% and 77% of the flux for the B, 
V and I band respectively regardless of the size and shape 
parameters of the sources. These ratios are used as the aper- 
ture corrections for the magnitudes of the observed GCs (see 

MM- 

Figure 13] shows a histogram of the FWHMs of the real 
GCs (secure sample) compared with the FWHMs of the sim- 
ulated objects (King profiles convolved with the instrument 
PSF). The excellent agreement between the two histograms, 
both at the "peak" and the "tail" of the distributions (the 
small discrepancies are most probably related to the fact 
that only one PSF model has been used for the whole mo- 
saic), confirms that the simulated objects are representative 
of a population of GCs spanning the whole range of shape 
parameters (i.e., Vc and c) expected for real sources. More- 
over, it ensures that the "tail" of the FWHM distribution 
is not due to contamination by spurious sources, but it is 
produced by real GCs (and precisely by the most extended 
and faint). The "tail" contains a significant fraction of the 
sources; for example, the fraction of GCs with FWHM be- 
tween 3 and 4 pixels equals 10% of the whole sample. In 
i]3.1l we reported that we imposed a FWHM limit of 4 pixels 
when selecting GCs. Based on the results of the simulation, 
we remark that choosing a more strict limit, although help- 
ing in excluding non-GC objects (e.g. galaxies), would also 
significantly reduce the number of real GCs in the candidate 
list. 

We further tested the FWHM limits by simulating the 
most extreme faint and extended GCs and measuring their 
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Figure 2. Results of tlie artificial source simulation for the V (top) and I (bottom) bands. The results were screened using the same 
criteria used to define the "GC candidates" sample ( ^3.11 1. Left.- Detection rate (i.e. percentage of detected sources over simulated sources) 
as a function of the input magnitude. The shadowed areas delimit the 90% and 75% completeness regions. The bold line represents the 
detection rate for the simulation run between the detection limit {Rqc = 25 ") and the galaxy R25- The light gray lines represent the 
results from simulations run within limited galactocentric ranges. The dependence of the detection rate curves to the radial distance from 
the galaxy center has been properly taken into account when correcting the V magnitude histogram for the effects of incompleteness 
(see i|4.1l l. Right - Ratio of the aperture flux measured by SExtractor over the total input flux of the simulated source, as a function 
of the input magnitude. The median value for the data points over 75% completeness, represented by the dashed line, has been used to 
compute the aperture corrections applied to the fluxes of the real GCs (see i|2.2l l. 



size. We defined the limiting magnitude for this simulation 
as that corresponding to the 90% completeness (calculated 
with the method described above, but excluding the FWHM 
criterion). We set the c parameter to 2.0 and we chose a rc 
of 4 pc (0.26 pixe ls) , which is in the u pper range of observed 
core radii (see [Jordan et al.lbOOSi '). The limiting FWHM 
found with this method was ~4 pixels, in agreement with 
the adopted selection criteria. 



4 PROPERTIES OF THE GLOBULAR 
CLUSTERS 

4.1 V band Luminosity Function 

The GC populations of almost all galaxies show a remark- 
able characteristic: their Luminosity Function (GCLF), plot- 
ted in magnitude units, is commonly represented by a Gaus- 
sian peaking at My ~ niag, with a dispersio n of 
av ~ 1.2 mag (e.g. lHarrij[l99ll : lAshman fc Zepj|l998l ). al- 



though minor deviations from this shape may be expected 
for a number of reasons (e.g. [Jordan et al.l |2007| . Villegas 
2010 for an extensive discussion on Gaussian representation 
of the GCLF in early type galaxies). A formation scenario 
that would justify a Gaussian distribution does not exist. 
Indeed, the GCLF has been tested agai nst other typ es of 
distributions, e.g., the student function l|Seckeilll992l '). For 
this reason, and given the skewness of the GCLF, the peak 
of the distribution is usually referred to as "turnover". The 
location of the peak appears to be fixed and therefore, when 
expressed in terms of the appare nt magnitude, it can be use d 
to estimate the galaxy distance (|Kundu fc Whitmorell200ll ). 

The left panel of Figure [4] shows the V band histogram 
for the candidate GCs (light gray), the secure sample of GCs 
(dashed shadow) and contaminating stars (gray shadow). 

The expected contamination due to foreground stars 
has been evaluated using the TRILEGAL online toojfl 

^ http:/ /stcv. oapd.inaf.it/cgi-bin/trilegal 
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Figure 3. Comparison of the FWHMs (in units of WFPC2 pix- 
els) of the simulated sources (solid line) versus the FWHMs of the 
measured sources in the data (shaded histogram) for the V band. 
The histograms are normalized to the total number of sources 
with FWHM between and 4 pixels. This result is valid for the 
"GC candidates" sample (no colour selection was applied in this 
simulation). 

l|Girardi et al.l [20051 ). which simulates the spatial distribu- 
tion of stars in the Milky Way. Assuming models for the 
star formation rate, age-metallicity relation and initial mass 
function of the Galaxy, TRILEGAL computes the age, metal- 
licity, mass and apparent photometry for the expected pop- 
ulation of Galactic stars along the line of sight within the 
desired field of view. We queried the tool using the default 
parameters (Chabrier function for initial mass function, ex- 
ponential disk model for extintion, thin disk -I- halo -I- bulge 
components for galaxy structure), for HST/WFPC2 (Vega 
system) magnitudes down to m,^ — 28 mag and for a field 
equivalent to the size of the D25 of NGC4261. Inspecting 
the histogram (Figure^)) we can see that the foreground star 
contamination is negligible. Moreover, this estimate should 
be considered as a rigorous upper limit, since this popula- 
tion was not filtered according to the selection procedure 
applied to the GCs ( ^3TT|) . 

4-1.1 Fitting The Luminosity Function Turnover 

We fitted a Gaussian profile to the histogram of the secure 
(colour selected) GC sample (shown in Figure |4]) in order to 
derive the peak location and FWHM of the GCLF. 

Notice that the distributions appear as skewed Gaus- 
sians due to the cutoff imposed by the incompleteness at 
faint fluxes. In particular, due to the colour selection and 
the more significant effect of the incompleteness in the / 
band data (caused by the higher background; see Figure [2]), 
the faint end of the GCLF of the secure GC sample might be 
biased by the high / band detection threshold. If this were to 
be the case, the result would appear as a deficit of blue ob- 



jects at faint magnitudes. Instead, comparing the colour his- 
tograms of objects brighter and fainter than >24 mag, 
we find that such a selection effect is minimal. 

Nevertheless, in order to obtain a conservative result, 
we estimated the magnitude at which such incompleteness 
starts affecting the sample using the following argument. 
The faintest object detected in the / band, at the 75% com- 
pleteness level, has a magnitude ~ 24.0 mag. Therefore, 
according to the minimum colour of the blue GC subpopula- 
tion {V — I ~0.6 mag; see jjS.lf) . a faint blue GC is expected 
to have a V magnitude of -^24. 6 mag (which is close to the 
75% completeness level for the V band). We adopted this 
value as the limit at which we "truncated" the fit of the 
GCLF. 

The fit was performed using the Sherpa package v4.2.1 
(which is part of the CIAO tool suite v4.2). The number of 
sources per bin of the Luminosity Function (LF) was high 
enough to allow the use of statistics. The uncertainties for 
each bin were estimated assuming a Poissonian distribution. 
The errors we report refer to the 1-a (68.3%) confidence 
bounds from the Ax^ projection for 1 interesting parameter. 
The number of foreground stars and other contaminating 
objects has been accounted for by including a constant in 
the fit. The resulting value is compatible with the estimated 
number of contaminating stars per bin (see Figure |4]). 

The turnover of the GCLF is located at m^, — 
25.llog mag, and the Gaussian distribution has a FWHM of 
3.llo 8 mag (oi' = l-3lo 3 rnS'S)- Usin g the Third Reference 
Catalogue of Bright Galaxies (RC3; de Vaucouleurs et al.l 
[199 9) we verified that reddening effects on peak location 
are negligible {Ab-v = 0.02 mag). The large errors on the 
peak location and FWHM are a consequence of the neces- 
sity of truncating the fit around the turnover. Neverthe- 
less, the best fit FWHM value is in close agreement with 
the results for the sample of elliptical galaxies studied by 
iKundu fc Whitmord (|200lh using data obtained with a sim- 
ilar configuration ( WFPC2, filters F555W and F814W). The 
results of the fits are shown in Figure 3] and Table [2] 

The turnover of the GC LF has been proven to b e 
a reliable distance indicator l|Kundu fc Whitmord I2OOII ). 
Adopting a A'/{?'^ = -7.4 mag for the peak of the GCLF 
(with an error negligible in comparison to our photo- 
metric errors) we derived a distance modulus m — M = 
32.5 

io.6 mag- The corresponding dista nce is 31.6lsy Mpc , 
in agreement with the estimate by IJensen et al.l (|2003l ) 
(29.4 ± 2.6 Mpc) bas ed on surface brightness fiuctuations. 
ITuUv fc Fisher! (Il988l ) measured a heliocentric velocity v 
= 2202 ± 75 km/s for NGC 4261, accounting for the in- 
fall of the Local Group to the Virgo cluster according to 
the model described in ITuUv fc Shaval (|l984l ). This veloc- 
ity corresponds to a distance of 31.3 ± 1.2 Mpc (adopting 
an = 70.4 ± 1.4 km/ s/Mpc - Seven Years WMAP Ob- 
servations; |jarosi^^eF|alJi2011.), which is consistent with our 
measurement at the 1-cr level. 

4- 1-2 GC Specifie Frequency 

The GC Specific Frequency (SFs; number of GCs per unit lu- 
minosity) of a galaxy is defined as Sn = Nac x ioO-i(Mv+i5) 
(|Harris fc van den Ber"ghlll98ll ). 

An estimate of the total number of GCs (Ngc) within 
the D25 was determined for both the secure and the candi- 
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Figure 4. V magnitude histograms for the GCs candidates and Gaussian fit. Left. {V magnitude histograms) - Candidate GCs are 
shown in hght gray, the subsample of candidate GCs selected in colour (secure GCs) with a dashed shadow and the contaminating stars 
(from the th.ilegal simulation) in gray shadow. Right. (Incompleteness-corrected V magnitude histogram and fits) - Candidate GCs 
are shown in light gray, the subsample of candidate GCs selected in colour (secure GCs) with a dashed shadow. The light gray area 
represents the distribution of the candidate GCs after applying the incompleteness correction, performed as described in i|3.2l The fit 
refers to the secure sample (green dashed line) and to the corrected candidate sa mple (violet dashed line) . The bin sizes have been 
optimized according to the Freedman-Diaconis criterion jPreedman fc DiaconisiriOSlI) . using the mean photometric error as the minimum 
allowed value. 



Table 2. Gclf Fit Results 



Sample 


Peak 
[mag] 


FWHM 

[mag] 


Amplitude^ 


Constant^ 


Integrated Number 
of Sources 


Distance 

[Mpc] 


Sjv 


Secure GCs 
Candidate GCs 


25 1 + 1-0 
^^■1-0.6 
25 5+0-2 


o 1+0.6 
•^■'^-0.8 

'^•'-'-0.3 




1-^-1.2 

t: 9+3.2 


1242 
2363 


31.6tJ« 
38.0I35 


1.5±0.3 
2.8±0.5 



t Value of [central] bin. Bin size: AV = 0.27 mag. 



date GCs by fitting tiieir histograms witii a Gaussian and 
a constant -which accounted for contaminating objects. Tlie 
LF for the candidate GCs sample had been previously cor- 
rected for incompleteness effects by splitting the sample into 
the three galactocentric ranges defined in Figure (2] applying 
the corresponding detection rate curve, and then stacking 
back the results into one histogram, shown in Figure 21 The 
parameters of the best fit Gaussian are shown in Table (2] 

Notice that the distance modulus derived from the fit 
of the LF for the candidate GCs (see Tableland is 
evidently overestimated. The peak of the Gaussian is driven 
towards fainter fluxes where the contamination by spuri- 
ous sources is amplified by the high completeness correction 
factor applied. Unfortunately, these sources cannot be dis- 
tinguished from real GCs using the FWHM and elongation 
criteria. For this reason, the total number of GCs resulting 
from the candidate sample should be considered as un up- 
per limit. On the other hand, the number derived from the 
secure sample, which was not corrected for incompleteness, 
represents a lower limit 

The integral of the GCLF models yielded Nqc ~1242 
and Nqc ~2363, for the secure and candidate sample re- 
spectively. Assuming an integrated magnitude for the galaxy 



of = 10.01 ± 0.06 mag (deriv ed as described be low) 
and a distance of 29.4 ± 2.6 Mpc (| Jensen et all 120031 ') . we 
calculated that: Sn ~1-5 ± 0.3 for the secure sample and 
Sjv =2.8 ± 0.5 for the candidate sample. Both numbers are 
within the range of typica. 1 values for elliptical galaxies (e.g. 
iKundu fc Whitmorell200ll ). 

The galaxy asymptotic magnitude was determined from 
a 2 D fit performed on the V band WFPC2 data using GAL- 
FIT (|Peng et al.ll2002h . This software is able to fit the surface 
brightness distribution of a galaxy to a 2D model account- 
ing for the effect of PSF blurring. We chose a Sersic model 
(plus a constant component accounting for the sky fiux) and 
used the SEXTRACTOR results for the galaxy as initial fit 
parameters. We used a Tiny Tim PSF image (generated as 
described in ij3.2|) . while we let GALFIT produce the weight 
image independently. We also applied a mask to exclude the 
central dust ring from the fit. The best-fit model has a Ser- 
sic index of 5.7, effective radius of ~70", major axis Position 
Angle (P.A.) of -21 deg and axis ratio (b/a) of 0.8. The ra- 
tio of the sum of residuals over the model-integrated light 
is less than 2% for all the filters, indicating that the model 
is a good representation of the starlight distribution of this 
galaxy. The galaxy magnitude was derived from the integral 
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counts under th e model applying the conversion formulae of 
iHoltzman et alj (|l995ah . 

The model-subtracted image showed a residual light 
ring around the galactocentric radius Rgc ~ 15", while we 
did no t find evidence for the tidal arm reported bv lTal et alJ 
l|2009 h (probably due to the lower exposure time). The pecu- 
liar pattern of the residuals is most likely related to the boxy 
nature of the galaxy. GALFIT is able to fit a boxy model with 
a fixed boxiness, while probably the isophotes of NGC 4261 
have a radial variation in this parameter (such behavior 
has b een reported in many ellipticals; e.g. iNieto fc Bender! 
I1989D . 

4.2 Colour Distribution 

An important result regarding GCs in massive galaxies is 
their bimodal col our distribution (e.g. lAshman fc Zepj|l99^ : 
iPeng et al.ll200d ). generally attributed to a metallicity dif- 
feren ce between two subpopulations (e.g. iBrodie fc Stradeil 
I2OO6I ). There are different scenarios which can account for 
such bimodality. The best supported is probably the merg- 
ing scenario, according to which the bluer (i.e., metal poorer) 
GCs are "donated" by merging galaxies, while the red (i.e. 
metal richer) GCs are formed in the merger, along wit h 
the bulk of the galaxy field stars (| Ashman fc Zepj|l992l ). 
This model has the advantage of naturally explaining the 
subpopulations and their different radial distributions (see 
Section iJS}, but does not address the extent to which the 
newly formed GCs can account for the significantly higher 
SF measured in ellipti cal galaxies w ith respect to those of 
the spiral galaxies fe.g. lHarrislll99ll ). which are supposed to 
be their progenitors. It is also not yet clear whether a signif- 
icant fraction of the young star clusters formed in merging 
galaxies can survive to account for the red GCs population 
of ellipticals. Valid alte rnatives are the multi-phase diss i- 
pational collapse model l|Forbes. Brodie. fc Grillmaiill 19971 ). 
which assumes that blue and red GCs are formed at different 
stages of the galaxy formatio n, and the accretion scenario 
l|Cote. Marzke. fc WestjllOgsl ). in which the red GCs repre- 
sent the intrinsic GC population of the galaxy while the blue 
GCs are acquired from lower mass galaxie s during dry merg- 
ers or by tidal stripping. As reviewed by iBrodie fc Stradeil 
l|2006r) . the three models are not mutually exclusive, and 
their effects may sum to generate the observed GC subpop- 
ulations in ellipticals. 

Figure [5] shows the V-I histogram for the secure GCs 
sample. The GC colour distribution doesn't show the clear 
bimodality typic al of the GC colour dis tribution in early- 
type galaxies (e.g. IBrodie fc Straderl2006l ). In order to quan- 
titatively discern the number of components in the distri- 
bution, we perfor med a mixture model analy sis using the 
NMIX method of iRichardson fc GreenI l|l997h . which im- 
plements a Bayesian test of univariate normal mixtures. 
We run the code f o llowin g the suggestions reported in 
IRich ardson fc GreenI (Il997l ) about the priors and hyperpri- 
ors (namely, choosing ^ = 1.1, K = 2, h = 2, and setting the 
rest as default), even though the test was demonstrated to 
be robust with respect to the choice of the initial parame- 
ters. The test did not show significant evidence for a bimodal 
GC colour (Bayes factor of 1 and 0.5, for the unimodal and 
bimodal models respectively). However, the range of colours 
of our GC sample spans the typical colours of both the blue 



Table 3. V-I Fit Results 



DOUBLE COMPONENT 


Component 


Amplitude^ 


Peak 
[mag] 


FWHM 

[mag] 


Blue Gaussian 
Red Gaussian 
Constant 




, 97+0.06 
^•^'-0.08 


30+°-^ 


Xl = 1.11, 10 d.o.f. (P(Q_„„,„e) = 35%) 


SINGLE COMPONENT 


Component 


Amplitude^ 


Peak 
[mag] 


FWHM 

[mag] 


Gaussian 






40+0.02 



Xi = 1.06, 14 d.o.f. (P(Q_„,,„,) = 39%) 



t Value of [central] bin. Bin size: AV — I = 0.06 mag. 

and red GCs in elliptical galaxies, while it is too wide to arise 
from a single GC population, suggesting that most probably 
the two subpopul ations are b l ended in our data. Moreover, 
as pointed out bv lPeng et al.l (|2006l ) in their work on early- 
type galaxies in the Virgo cluster, even for those galaxies in 
which the relatively low number of GCs does not allow one 
to distinguish between unimodal and bimodal colour distri- 
butions, the result of a decomposition into two components 
provides a reliable description of the system. For these rea- 
sons, we assume that the distribution consists of superim- 
posed blue and red subpopulations and so we fit the colour 
histogram using two Gaussian components, plus a constant 
to account for the contamination by foreground stars (esti- 
mated as 1-2 objects per bin using the TRILEGAL results 
- see M.l^ and background galaxies. The uncertainties for 
each bin were estimated assuming a Poissonian distribution. 

Table reports the result of the fit. For completeness, 
the table also shows the result from a single Gaussian fit. In 
the double-component case, the large errors in the Gaussian 
normalizations refiect the fact that the two amplitudes are 
strongly correlated. Nevertheless, the locations of the peaks 
are accurate enough to allow us to identify a blue and a red 
subpopulation, with average colours T/-/=1.0ll° °^ mag and 
1/-/=1.27+Q°g rnag respectively. These value s are in good 
agreem ent (within Icr) with those reported bv lLarsen et al.l 
(l200ll ) based on a similar analysis of the GC population of 
a sample of nearby early-type galaxies, observed using deep 
exposures in the F555W (V") and F814W (/) HST fihers. 

The best-fit model along with the individual compo- 
nents are shown in the top panel of Figure \S\ The bottom 
panel shows the residuals with respect to the best-fit model 
in units of standard deviation. 

We arbitrarly divided the secure GCs into blue and red 
by setting the threshold at V-I — 1.15. The discriminating 
value is representative of the 50% probability for a GC to 
belong to the blue or red subpopulation (the intersection of 
the Gaussians in Figure [SJ. This value is consistent with the 
separation of the two subpopulations in other galaxies (e.g. 
I Ashman fc Zepilll998l '). 



© 0000 RAS, MNRAS 000, 000-000 



Studying the asymmetry of the GC population of NGC 4261 9 



110 
100 

90 
80 

i ^° 

1 60 
' 50 
! 40 
■ 30 
20 
10 


S 2 



1 secure GCs 
double gaussian + const 

■ blue component 

■ red component 




0.5 



1.0 



1.5 



V-l [mag] 



Figure 5. V-I colour histogram for the secure GCs sample and 
fit. The data have been fit using two Gaussians plus a small 
constant used to account for contamination. In the figure, we 
show the two Gaussian components (blue and red dotted lines) 
and the total convolved model (gray dashed line). The addi- 
tive constant is not showed in the figure. The bin sizes have 
been optimized according t o the Freedman-Diaconis criterion 
dFreedman &: Diaconij[l98ll) . using the mean photometric error 
as the minimum allowed value. 



In order to estimate the relative abundance of the blue 
and red subpopulations, we performed the following Monte- 
Carlo simulation. We sampled each of the parameters of the 
two-Gaussian fit (described above). The sampling distribu- 
tion was a multi-variate Gaussian with a mean equal to the 
best-fit value of each parameter and a FWHM derived from 
the covariance matrix of the fit. At each trial, the set of sam- 
pled parameters defined two Gaussians whose flux fractions 
with respect to the total flux from the combined model were 
recorded. Finally, the expectation value for the fraction of 
the blue/red component was calculated as the median over 
all the trials. With this method, we estimated that the blue 
(red) GCs represent the 60^2^% (40^2^%) of the total sam- 
ple of 718 GCs. The upper and lower uncertainties on the 
fraction represent the 68.3% quantile (1 cr) around the me- 
dians. 



5 SPATIAL DISTRIBUTION OF SOURCES 



Prompted by the fi ndings of IZezas et al.l (|2003l ) and 
iGiordano etal] (|2005l ). we investigated the spatial distribu- 
tion of the secure GC sample. 

The region of the galaxy within Rgc < 25" was ex- 
cluded from the analysis due to the high background which 
drastically limited the source detection sensitivity (see i]3.2|) . 
Since the GC detection rate depends on the galactocentric 
distance (Figure [2)| and affects the relative number of de- 
tected sources at each radius, introducing artificial asymme- 
tries, we decided to limit the analysis to those GCs having 
a magnitude brighter than the 75% completeness limit for 
the V band (corresponding to my = 24.6 mag). 



5.1 Radial Distribution 

In the top panel of Figure [S] we present the radial profiles of 
the blue and red GCs. As a reference, we also plot a scaled 
l/-band light profile model of NGC 4261 (eis derived from 
the GALFIT fit; see >i4.1.2p . The data points represent the 
surface densities (number of sources per arcsec'^) over ellip- 
tical annulii centered on the galaxy center. The orientation 
and ellipticity of the annulii have been chosen to match the 
P.A. and ellipticity of the galaxy light profile. The profiles 
have been normalized to match at the first point. 

We used funtools (jMandel. Murray, fc Rolil bOOll ) to 
count the number of sources within each annulus. The an- 
nulii sizes (i.e. binning of the plot) have been chosen in order 
to contain a statistically significant number of sources (at 
least 20 sources per bin) that would allow the use of a 
test. The uncertainty for each bin is simply the square root 
of the counts divided by the annulus area. 

The distribution of the overall GC population looks con- 
sistent with the trend of the light profile, at least from the 
galaxy center up to R/R25 = 0.75, where it starts declin- 
ing, mostly driven by the drop of the red GC component. 
The difference between the trends of the blue and the red 
GCs radial profiles is more evident when the distributions 
are plotted in terms of Specific Frequency (Sn), calculated 

as SN{i) = N'gc X 10" ''("^+^5' , where i is the annulus 
index (bottom panel of Figure [SJ. In fact, recalling that Sn 
expresses the relation between the number of clusters and 
the galaxy light, the systematic decrease of S^^" denotes a 
clear inconsistency with the galaxy surface brightness. On 



the other hand, 5*. 



varies around a constant value. 



Such behaviour seems to be in contrast with the com- 
monly accepted scenario (see according to which the 
red (young) GCs follow the galaxy l ight, while the blue (old 
are spread almost un iformly (e.g. iBrodie fc Straderl 1200' 
lAshman fc Zepj|l998l ). 

On the other hand, the strong blending of the colours of 
the two populations of our sample (see may not have 

allowed us to separate properly the two profiles. 

A Kolmogorov-Smirnov (K-S) test has been used to 
compare the radial distribution of blue and red GCs against 
a Sersic brightness profile (whose parameters have been de- 
rived from the GALFIT fit; see ^4.1.2(1 . Since the Sersic pro- 
file is defined assuming circular symmetry, the dependence 
of the galactocentric distance of a GC {Rgc) on the ellip- 
ticity of the galaxy had to be removed as well. The variable 
used to perform the test was the major axis a of the ellipse 
intercepting the GC and defined by the same center, elon- 
gation and P.A. as for the galaxy. This can be calculated 



RGC,i X v/e2 cos2(0, - P.A.) + sin2(e(, - P.A. 



(5) 



where e is the elongation, 9 the azimuthal angle {6 — 
towards North and increases counterclockwise) and i the GC 
index. The test returned a probability well below the 0.1% 
level for the blue GCs and a probability of 0.3% for the red 
GCs, therefore rejecting the null hypothesis that their radial 
distributions follow the galaxy fight profile. 
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Figure 6. GCs radial profiles. Top. Comparison of the source 
densities of the GCs sample (defined in Section Jsjl and of the blue 
and red subpopulations (green, blue and red lines respectively) 
with the scaled light profile model of the galaxy (violet line). The 
data points represent the source densities evaluated over elliptical 
annulii. The radial coordinate is expressed in terms of the galaxy 
semi-major axis R25 = 2.04". The vertical error bars represent 
Poissonian errors. The profiles have been normalized in order to 
match at the radius of the first annulus. Bottom. Radial profile of 
the Specific Frequency (Sjv) of the blue and red subpopulations 
(blue and red lines respectively), calculated as S'jv(j) = Nq^-, x 
-|^q0.04(m^+15)^ where i is the annulus index. 



5.2 Azimuthal Distribution 

We evaluated the azimuthal distribution of the overall GC 
sample as well as of the blue and red subpopulations by 
comparing their source densities against the galaxy surface 
brightness within a set of elliptical wedges centered on the 
galaxy nucleus. The parameters of the ellipses (center, ellip- 
ticity) have been obtained from the GALFIT fit (see H4.1.2p . 
The angular extent of the wedges has been chosen to ensure 
Gaussian statistics for the majority of the bins (>20 sources 
per bin). 

In our mosaics, 4 areas within the galaxy D25 were cov- 
ered by the WFPC2 PC camera (field of view of 35"x35"). 
We had to exclude these locations from the current analy- 
sis due to their different sensitivity (which would infiuence 
the surface brightness estimation) and higher noise (which 
could create artificial deficits of GCs). When measuring the 



galaxy fight, we also masked out all the GCs (secure sam- 
ple), foreground stars and background galaxies. 

The measurements were repeated for different rotations 
(offsets) of the ellipse wedges in order to address the signif- 
icance of any density excess. The wedge offset is defined on 
the rotation of a complete set of wedges measured counter- 
clockwise with respect to the North. For each offset, a test 
between the azimuthal profiles of the GC density and that 
of the galaxy fight has been performed in order to identify 
the wedge rotation for which the anisotropy is maximized. 
Notice that, in this context, the azimuthal surface bright- 
ness distribution of the galaxy simply represents a uniform 
comparison distribution which also includes the modulation 
due to the geometry of the wedges alone. 

Figure [7] shows the source density azimuthal profile of 
the whole GC sample as well as the blue and red subpopu- 
lations for the wedge-offset that maximizes the statistic 
for each series. The solid lines represents the best-fit of the 
galaxy light to each data series (used to calculate the x'^)- 
The counter-intuitive modulation of the azimuthal distri- 
bution of the galaxy fight is due to a combination of two 
effects: the boxiness of the galaxy (so that the wedges along 
the major axis have larger areas but include less flux), and 
the masking we applied. The bottom panel of the figure re- 
ports the fit residuals expressed in terms of a with respect 
to the fitted model (i.e. starlight distribution). 

The plot shows a clear density enhancement (more than 
2 a) at P.A. ~225±18°, and significant deficits at P.A. 
~120±18° and at P.A. -270±18°. The peak-to-peak varia- 
tion is of the order of several a\ the red GC subpopulation 
seems to follow the same modulation, while the statistics for 
the blue distribution are too poor to derive any conclusion. 

The GCs in the overdensity regions do not show any 
difference in brightness or colour distribution with respect 
to the "field" GCs (the KS test can not distinguish the two 
distributions), their V and V — I histograms closely resem- 
bling the ones of the whole secure GCs sample (shown in 
Figures [5] and lU . 

This test showed that the reduced (x^) for the whole 
sample was as high as 2.5 for 9 degrees of freedom. Since we 
deliberately select the highest the usual x? probabil- 
ity density function will not provide a correct confidence 
level. Instead, we had to compare this value against the x? 
generated by a population of GCs distributed according to 
the galaxy light, in order to exclude the possibility that the 
overdensities of GCs are statistical fluctuations. To repro- 
duce the proper xS probability distribution, we performed 
the following test. We generated n points (where n equals 
the number of wedges) distributed according to a Poisso- 
nian density function around the average number of GCs 
found within a wedge. Such a set of data represented the 
equivalent of the GCs azimuthal proflle measured for a cer- 
tain wedge rotation. The data were then fit with a constant. 
The correct way to run this test would have been to use the 
galaxy light azimuthal distribution to derive the average of 
the Poissonian distribution in each wedge, but, since we as- 
sume a constant Nqc to light ratio in each wedge, the two 
methods are exactly equivalent (in the latter case the ap- 
propriate fit would just be a rigid shift of that curve along 
the y-axis). We simulated a set for each wedge offset and we 
recorded the maximum x^ among the sets, as for the real 
data. 
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Figure 7. Source density azimutlial profile. Tlie data points 
represent tlie density of GCs witliin wedges of 36° for the GCs 
sample (green diamonds) and for the red and blue subpopula- 
tion (red dots and blue squares respectively) as a function of the 
wedge P.A., for the wedges offset maximizing the anisotropy of 
each series (P.A. = towards North and increases counterclock- 
wise). The solid lines represent the best fit of the galaxy light 
to each data series. The bottom panel of the figure reports the 
scatter of the three data series with respect to the fit in units of 
<T. The data points have been slightly shifted along the x — axis 
for presentation purposes. 



The whole simulation was repeated 2000 times in order 
to generate the xl probability distribution. The histogram 
in Figure [8] shows the resulting distribution of the simulated 
max xi- The arrow shows the value of the xt measured 
from the data. The integral of the simulated xi exceeding 



this value corresponds to a probability P^2^^2 < 0.1%. 
This result indicates that the asymmetry in the azimuthal 
distribution of the GCs with respect to the galaxy hght is 
statistically significant at the 99.9% level. 

As a further test, we compared the distributions of 
the candidate GCs, and the blue and red subpopulations, 
against a uniform distribution using a K-S test. The KS 
test results did not allow us to reject the null hypothesis (i.e., 
that the azimuthal profile follows a uniform distribution), for 
either of the three GCs subpopulations (we obtained prob- 
abilities of 7%, 35% and 73% for the whole sample and the 
blue and red subpopulations, respectively). In order to assess 
the significance of this result, we calibrated the K-S test in 
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Figure 8. Probability distribution of xi produced by the fit of 
a constant to random sets of 10 data points distributed according 
to a Poisson distribution. The simulation is meant to reproduce 
the xi derived for a uniform azimuthal distribution (see text for 
more details). The arrow represents the xi obtained fitting the 
azimuthal distribution of the GCs to the galaxy light. The area 
covered by xi equal or higher than this value corresponds to 
a probability P 2~,^2 < 0.1%, therefore indicating that the 

asymmetry in the azimuthal distribution of the GCs with respect 
to the galaxy light is statistically significant. 

a similar fashion as for the x^ method. For each population 
(secure GCs, red subpopulation and blue subpopulation) , 
we generated m random angles (where m is the number 
of sources in the population) uniformly distributed in the 
range [0, 2n]. The cumulative distribution of this sample was 
compared against the cumulative of a uniform distribution 
and the maximum distance D between the two was recorded 
(such distance is in fact the K-S statistic). The simulation 
was repeated 10000 times in order to generate the probabil- 
ity distribution for D. The integral of the values exceeding 
the D measured in the K-S test on the real data correspond 



to the probabilities PoyD^^, ~ 10%, PoyD^^, ~ 70% and 
Pd->d ^ ^ 30% for the secure GCs sample and the red and 
blue subpopulations, respectively. This result confirmed that 
the K-S test is not sensitive enough to identify the differences 
seen between the samples. 



5.3 The two-Dimensional Distribution of the GC 
Population 

In order to further investigate any local enhancements in 
the spatial distribution of the GCs we used the CIAcQ tools 
csmooth and vtpdetect. 

These tools were originally designed to identify clusters 
of photons in X-ray images. For this reason, we created "im- 
ages" of the GC locations to be used by the tasks, with pixel 
value equal to 1 where a GC was present and elsewhere. 

The csmooth tool smooths the image using an adaptive 
Gaussian kernel: the scale for each area is increased until the 
counts under the smoothing kernel exceed the significance 
threshold over the (local) background set by the user. We 

■* The Chandra data analysis package, supported by the Smith- 
sonian Astrophysical Observatory. 
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Figure 9. Local enhancements of GCs. Left. Over a mosaic of NGC 4261 (/ band), we overlay the sectors over which we evaluated 
the GC azimuthal densities shown in Figure [7] (black ellipse). The inner and outer radius of the ellipse correspond to Rqc = 25" and 
Rgc = ^25 respectively. The elliptical annular wedges we show are rotated to the angles at which the anisotropics in the GC azimuthal 
distribution are maximized. The filled wedges map the density enhancements shown in Figure[7] the color intensity is proportional to the 
significance of the excess (i.e. value of the corresponding - positive - residual in bottom panel of Figure [T]!. Right. Overlaid on a mosaic 
of NGC 4261 (/ band), the shaded ellipses show local enhancements in the distribution of GCs detected with a Voronoi Tesselation and 
Percolation (VTP) algorithm. Each clustering is represented by an ellipse whose axes size equals the a of the spatial distribution of its 
GCs. The black shadows represent the areas of higher clustering of the secure GC sample as detected by the ciao task csmooth. The 
results from the VTP algorithm and from csmooth are clearly consistent with those from the wedges test, and they point to a GC density 
excess along a NE-SW direction. 



set a minimal significance threshold of 3 a. The background, 
which in our case represents the uniform component of the 
GCs distribution, has been estimated locally. The task is 
robust with respect to the modification of the other param- 
eters (such as the initial guess of the smoothing scale), which 
mostly affected the computational time. 

The result of this procedure is shown in Figure [H] We 
can clearly identify two regions where the clustering of 
GCs appears to be significant (the regions appear slightly 
curved since we excluded sources near the center of the 
galaxy). These enhancements in the source density are lo- 
cated around P.A.~30° and P.A.~245° (NE and SW re- 
gions), in agreement with the locations of the peak in the 
azimuthal distribution of GCs (Figure [7]). 

The vtpdetect task exploits the Voronoi Tessellation and 
Percolation (VTP) technique to detect connected sets of 
photons above a threshold density. In our case, each GC 
played the role of a distinct photon. The task "cuts out" a 
triangular area around each source and calculates the source 
density (number/area). The background model is estimated 
by fitting the low-end of the cumulative density function 



(CDF). A density cutoff (used to distinguish areas contain- 
ing a source photon from background) is derived from the 
residual between the CDF and the background model. Con- 
tiguous areas over the density cutoff are then clustered to- 
gether using a percolation algorithm. A source built in this 
way is considered significant if it exceeds the desired false 
source limit due to background fiuctuations. A detailed de- 
scription of the VTP theory is presented in the CIAO Detect 
Reference Manual (December 20060- 

In order to use vtpdetect on non X-ray data we had to 
deduce the task parameters through a trial-and-error ap- 
proach. We estimated a GC "flux cutoff" comparing the 
source density in the outer regions of the image against the 
average area per source. The result was compatible with the 
default value for the relevant task parameter {maxbkgflux 
— 0.8). We adopted the default values also for the other 
parameters related to the background estimation {mintot- 
flux = 0.8, maxtotflux — 2.6, mincutoff = 1.2, maxcutoff 

^ ' http: / / cxc.harvard.edu/ciao/download/doc/detect_manual/| 
vtp_theory.html 
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— 3), after verifying that they did not significantly affect 
the detection resuhs. We set limit — 0.001 and scale = 1; 
these values have been chosen in order to limit the source 
deblending. 

The la VTP ellipses defining source clustering are 
shown in Figure [§1 The results from the VTP test confirm 
the excess detected by csmooth. 



6 DISCUSSION 

Using HST WFPC2 data in the B, F and / band (HST 
filters F450W, F606W and F814W, respectively) we charac- 
terized the GC population of the nearby elliptical galaxy 
NGC 4261. We studied the Globular Cluster Luminos- 
ity Function (GCLF), deriving a total number of ~2363 
GCs, corresponding to a Specific Frequency (SF; number 
of GCs per unit magnitude) of 2.8±0.5. From the loca- 
tion of the peak of the GCLF, we derived a dist ance of 
31.6t . l'''^^ Mpc (assurning a Mv ~ —7.4 mag - e.g. iHarrij 
ll99ll : lAshman k Zepj liggs'l. The colour distribution can be 
interpreted as the superposition of a red and a blue sub- 
population, with average colour V — I = LOl^o oe mag and 
1.271q'o8 niag respectively. 

Most importantly, we used different techniques to study 
the two-dimensional distribution of the GCs sample and its 
blue and red subpopulations. As discussed in i]5.2l and ^5.31 
we discovered that the GC population of NGC 4261 shows 
evidence for an asymmetry in the azimuthal distribution 
above the 99.9% confidence level. A similar effect is noticed 
in the red GC subpopulation, but on at lower significance. 
The result for the blue GCs is not statistically significant. 
The origin for such anisotropy may reside in the evolution 
of the galaxy in the recent past. Next, we discuss possible 
scenarios that address to this asymmetry. 

(i) Major merging event. The effects of a major merg- 
ing event represent the most obvious candidate to explain 
such a significant asymmetry. NGC 4261 does not show 
signs of recent strong interactions (shells, ripples, arcs, tidal 
tails or other distortions), therefore we exclude a major 
merger within the last 1 Gyr, since this is the timescale 
on which interaction fea t ures fade away (e.g. lQuinnlll984l : 
iHernauist fc Quinnll 19881 . Il989l ') . The difference between the 
starlight and GC distributions implies that, if the GC asym- 
metry is related to a major merger event (whether new GCs 
have been formed in situ, donated or just been displaced), 
the relaxation timescale for the GC population is larger than 
for stars. This is a complicated issue since GC systems ex- 
tend beyond the main stellar body of the galaxy, and their 
dynamics at large galactocentric radii are affected by the 
dark matter (resulting in higher velocity dispersions for the 
GCs than for the stars). To our knowledge, there is no re- 
ported evidence for such a difference in the relaxation time 
so far. 

A major merging event may have triggered episodes of 
star formation within remnant tails along the line of sight. 
In this case, the enhanced GC density along the NE-SW 
axis of NGC 4261 may just be due to a projection effect. 
Studies of galaxies at an advanced st ate of merging, such as 
NGC 7252, NGC 3256 or NGC 3921 (|Knierman et al .1120031 : 
ISchweizer et al.lll99 6l show indeed that tidal tails host sig- 
nificant populations of proto-GCs. Therefore a certain frac- 



tion of GCs form in highly asymmetric distributions. The 
question is what fraction of these proto-GCs can survive 
tidal destruction and whether they can maintain their origi- 
nal spatial asymmetry during the subsequent dynamical evo- 
lution of the host galaxy. 

In addition, there are cases suggesting incomplete violent 
relaxatio n which would naturally result in non uniform dis- 
tribution. ISchweizer et al.l (|l996l ) showed that the candidate 
GCs of the protoelliptical NGC 3921 may have experienced 
the same (incomplete) violent relajcation as did the average 
starlight. Nevertheless, their data (although only partially 
covering the galaxy extent) show indications of asymmetry 
both in the radial and in the azimuthal distribution of the 
GCs. This suggests that spatial asymmetries of GC systems 
may indeed survive the relaxation process. 

Numerical simulations suggest that massive boxy ellipti- 
cals (such as NGC 4261) are most probably the result of dis- 
sipa tionless (dry) mergers between ellipticals of similar size 
(e.g. iKhochfar &: Burkert|[2005l : lNaab. Khochfar. fc BurkertI 
l2006l ). Since this scenario does not involve gas, or all the gas 
is converted to stars prior to the collapse, the number of 
newly formed GCs along tidal tails would be negligible. 

Moreover, according to the merging scenario, the dis- 
placed GCs should mostly be younger (i.e. metal-rich). Since 
we observe asymmetry in both the blue and red subpopula- 
tions (although the result for the blue GCs is not statistically 
significant), we may conclude that the former hypothesis 
does not apply in the case of NGC 4261. 

Given the forementioned considerations, we conclude that 
a major merging event may be responsible for the observed 
asymmetry only if the relaxation timescale for the GC pop- 
ulation could be proven to be larger than for stars. 

(ii) Minor merging event. As mentioned above, the ab- 
sence of fine structure in NGC 4261 excludes the possi- 
bility of a recent (~ 1 Gyr) major merging, but it does 
not rule out this possibility for minor mergers, since these 
events have weaker effects on the galaxy light distribu- 
tion. This is the case in the elliptical galaxy NGC 1052, 
which underwent a recent (1 Gyr) mergin g event with 
a gas rich dwarf (Ivan Gorkom etal] 1 19861 '). The galaxy 
shows a weak rotation about its minor-axis (suggesting tri- 
axiality) and significant minor-axi s pos ition angle swings 
(see notes iniDa vies fc Birkinshawlll988l ). but, similarly to 
NG C 4261, it exhibits a ver y low fine structure parame- 
ter ISchweizer fc Seitzeij[l99^ ). iPierce et all (|2005l 'l though, 
found no young GCs associated with the recent merger event 
of NGC 1052. This may indicate that no significant increase 
in the GC population should be expected even in case of gas 
rich minor merging. 

Non-relaxed GC distributions can be easily explained by 
GCs donated in a recent minor merging event but their num- 
ber densities would be too small to generate statistically 
significant asymmetries over the galaxy light background. 
On the contrary, we can not rule out the possibility that 
the "indigenous" GCs have just been displaced during the 
merger. 

(iii) Galaxy interactions. Galaxy interactions are 
known to provoke displacements of the GCs systems. 
For example, the flyby encounter of NGC 1399 with 
the nearby NGC 1404 disturbed the velocity structure 
of the outer regions of the GC system of NGC 1399 
iNapolitano. Arnaboldi. fc Capacciolil |2002| ). or even gen- 
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erated a "tidal stream" of intergalactic GCs ()Bekki et alj 
[2003,, '). A smoothed map of the displa ced blue GCs candi- 
dates around NGC 1399 is reported in lBassino et alJ l|2006l . 
Figure 9). 

Moreover, simulations by IVesperini fc Weinberg! (|2000l ') 
showed that the perturbations on stellar systems of elliptical 
galaxies due to flyby encounters can last well past the clos- 
est passage of the perturber. Similarly, we can expect that 
the asymmetric distribution of the GC system of NGG 4261 
may be a remnant from the recent interaction with one of its 
close by companions. In order to identify the candidate per- 
turber, we calculated the R MS radial velo city of the mem- 
bers of the NGC4261 group (|Garcialll993 ) and, considering 
a time between 1 and 2 Gyr for the remnants of the inter- 
action to disappear, we estimated that the perturber should 
be found between 15' and 30' off NGC 4261. Unfortunately, 
even considering only the brightest galaxies {g < 14 mag) 
in this range, there are too many candidates to uniquely 
identify the possible perturber. 

In this sense, a close encounter similar to the one that 
happened between NGC 1399 and NGC 1404 may have 
shifted some of the GCs of NGC 4261 from the NW-SE 
"poles" towards the NW-SE edge-on plane, with the galaxy 
nucleus preventing us from observing the distribution at 
the innermost radii. The GC systems extend much fur- 
ther than the host galaxy light, therefore the displacement 
may have affected only the outermost GCs without affect- 
ing the stellar orbits and hence the starlight. This event 
may also justify the peculiar drop of the GC radial pro- 
files at large radii, shown in Figure [5] (although the effect 
may be due only to statistical fluctuations). If the GCs are 
indeed distributed continuously along a NW-SE axis, they 
would lie along a plane pe rpendicular to the rotation axis 
ijDavies fc Birkinshawlll9 8d). therefore the overdensity could 
be simply associated with the rotation of the galaxy. 

In addition, some GCs may have been stripped out of 
the interacting galaxy. However, since the candidate per- 
turbers are significantly smaller than NGC 4261 (therefore 
they have poorer GC systems), this effect can not explain 
alone the azimuthal variation of the GC density. In fact, its 
variation can reach up to ~ 50% (distance of the GC density 
from the fit of the galaxy surface brightness; Figure [7|), re- 
quiring a number of stripped GCs comparable to the system 
of NGC 4261. We find that there is no evidence for a dif- 
ference between the properties (colour, FWHM, axis ratio, 
etc.) of the objects in the overdensity regions and those in 
the "field" . This implies that the number of GCs subtracted 
from the interacting galaxy is relatively small or if the ac- 
quired GC population is numerically significant, then it has 
similar properties. 



The displacement of GCs under the effects of galaxy in- 
teraction or mergers and estimation of the relaxation times 
for GC "particles" embedded in mixed baryonic and dark 
matter environments is a subject that needs to be explored 
in order to understand the distribution of the GC popula- 
tions in the increasing number of galaxies which show non- 
uniform GC systems. 



7 CONCLUSIONS 

We performed an analysis of the characteristics and spatial 
distribution of the GC population of the nearby elliptical 
galaxy NGC 4261 using deep HST WFPC2 data in the B, 
V and / band (HST filters F450W, F606W and F814W, 
respectively). Our photometry of the GC system resulted 
in a sample of 718 secure GCs, complete down to m^^ = 
23.8 mag (90% completeness level). We were able to derive 
the Globular Cluster Luminosity Function (GCLF) down 
to the limiting magnitude V = 24.6 mag. The distribution 
peaks at my = 25.11^ ^ mag, consistent with a distance 
of 31.6j^g'''7^ Mpc (a ssuming a peak My ~ —7.4 mag - e.g. 
lHarris|[l99ll : lAshman & Zcpf 199& ). We studied the colour 
distribution and probed the bimodal beh aviour expected 
in ell iptical galaxies for GC systems (e.g. lAshman fc Zepll 
The V — I colour distribution is compatible with the 
superposition of a blue and a red component with average 
colour 1.011q q6 nrS'g ^^'^ l-27lo o8 '^''■g respectively. We cal- 
culated a Specific Frequency (number of GCs per unit mag- 
nitude) of 2.8±0.5, therefore relating the GC abundance to 
the stellar mass. We analyzed the two-dimensional distribu- 
tion of the GC sample, and we found that it shows evidence 
for an asymmetry in the azimuthal distribution above the 
99.9% confidence level. A similar effect is noticed in the red 
GC subpopulation alone, but at a lower level of statisti- 
cal significance. We discussed the origin of this anisotropy 
in context of galaxy interaction and evolution. Our main re- 
sults regarding the GC spatial distribution are the following: 

(i) The GC radial profile follows the galaxy surface 
brightness especially at galactocentric radii R/R25 < 0.6. 
The red subpopulation profile is systematically steeper than 
the starlight, while the profile of the blue subpopulation 
looks more consistent with it (Figure (Sjl. 

(ii) The azimuthal distribution of the GC population 
shows significant enhancements along a NE-SW direction 
(Figure [7] and [9]). The red subpopulation seems to follow the 
same modulation, reflecting each local enhancement /deficit. 
The result for the blue GCs is statistically not significant. 

(iii) We suggest that the peculiar asymmetric distribu- 
tion may be related to a past (dry) merger or interaction 
event. In particular, we favour the hypothesis that a fly-by 
encounter displaced the outermost GCs. The hypothesis of 
a past (older than 1 Gyr) major merging event cannot be 
ruled out, modulo the assumption that the GC relaxation 
time is longer than for the stars. We remark that the GCs 
may be distributed on an edge-on plane, since the galaxy 
nucleus prevents us from observing the distribution at the 
innermost radii. If this were to be the case, we cannot ex- 
clude that the overdensities could be rotation-related. 

The study of the anisotropy of the GC population of 
a galaxy is a relatively unexplored field, which offers the 
unique opportunity of exploiting the GCs to investigate the 
history of the host galaxy. Spectroscopic data are critical to 
investigate whether the GCs of NGC 4261 show any peculiar 
kinematics or they have typical dispersion velocities for ellip- 
tical galaxies. In particular, spectroscopy of the GCs in the 
overdensity regions will clarify whether they are displaced 
GCs and move as an independent system from the overall 
GC population. Finally, we enphatize the need for numeri- 
cal simulations aimed at modelling the evolution of the GC 



© 0000 RAS, MNRAS 000, 000-000 



Studying the asymmetry of the GC population of NGC 4261 15 



spatial distribution after a galaxy interaction or a merger 
and at estimating the relaxation times for GC "particles" . 
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Table 4. Catalogue Of Detected GCs 05 



GC 


RA (J2000) 


Dec (J2000) 


e 


FWHM 


Stellarity 


(S/N) 


B 


V 


I 


B-V 


V-I 




(ID) 


[hh:mm:ss] 


[dd:mm:ss] 




[pixel] 






[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


ta 


1 


12:19:24.991 


+05:47:39.48 


1.07 


1.78 


0.97 


278.1 


23.52+0.02 


22.67+0.06 


21.69+0.01 


0.94+0.01 


0.95+0.01 


mfin 


2 


12:19:22.278 


+05:47:40.68 


1.15 


1.68 


0.94 


32.4 


26.13+0.07 


25.28+0.06 


24.08+0.02 


0.94+0.05 


1.25+0.04 


3 


12:19:25.512 


+05:47:41.86 


1.06 


1.86 


0.98 


289.8 


23.54+0.02 


22.61+0.06 


21.52+0.01 


1.01+0.01 


1.10+0.01 




4 


12:19:24.360 


+05:47:42.57 


1.12 


2.15 


1.00 


89.5 


24.87+0.03 


24.13+0.06 


23.34+0.01 


0.84+0.02 


0.67+0.02 


5 


12:19:22.175 


+05:47:42.28 


1.18 


2.52 


0.97 


34.1 


25.61+0.04 


25.24+0.06 


24.33+0.03 


0.49+0.04 


0.84+0.04 




6 


12:19:26.772 


+05:47:44.50 


1.31 


2.05 


0.97 


42.4 


- 


25.04+0.06 


24.25+0.02 


- 


0.68+0.04 




7 


12:19:22.008 


+05:47:44.76 


1.22 


3.65 


0.97 


29.2 


- 


25.43+0.06 


24.44+0.03 


- 


0.96+0.05 




8 


12:19:24.579 


+05:47:45.71 


1.03 


1.91 


0.95 


103.5 


25.00+0.03 


23.95+0.06 


22.75+0.01 


1.12+0.02 


1.25+0.02 




9 


12:19:26.670 


+05:47:47.15 


1.11 


2.28 


0.98 


290.4 


23.56+0.02 


22.62+0.06 


21.54+0.01 


1.02+0.01 


1.08+0.01 




10 


12:19:22.306 


+05:47:46.92 


1.07 


1.83 


0.97 


243.8 


23.68+0.02 


22.82+0.06 


21.75+0.01 


0.95+0.01 


1.07+0.01 




11 


12:19:25.243 


+05:47:46.83 


1.05 


2.43 


0.98 


49.9 


25.33+0.04 


24.81+0.06 


23.91+0.02 


0.62+0.03 


0.84+0.03 




12 


12:19:25.279 


+05:47:50.12 


1.29 


1.93 


0.98 


557.4 


22.45+0.02 


21.58+0.06 


20.47+0.01 


0.96+0.01 


1.13+0.01 




13 


12:19:25.400 


+05:47:49.67 


1.04 


2.18 


0.99 


663.6 


22.19+0.02 


21.28+0.06 


19.98+0.01 


0.99+0.01 


1.39+0.01 




14 


12:19:27.281 


+05:47:49.17 


1.36 


3.03 


1.00 


123.3 


24.76+0.03 


23.77+0.06 


22.80+0.01 


1.08+0.02 


0.92+0.01 




15 


12:19:26.721 


+05:47:49.01 


1.01 


1.89 


0.95 


99.8 


24.94+0.03 


24.02+0.06 


22.89+0.01 


1.00+0.02 


1.16+0.02 




16 


12:19:24.447 


+05:47:51.17 


1.15 


2.22 


1.00 


158.7 


24.37+0.03 


23.41+0.06 


22.49+0.01 


1.05+0.01 


0.85+0.01 




17 


12:19:20.828 


+05:47:50.86 


1.41 


2.23 


0.95 


88.7 


25.31+0.04 


24.16+0.06 


22.98+0.01 


1.22+0.03 


1.22+0.02 




18 


12:19:21.749 


+05:47:52.19 


1.12 


1.89 


0.97 


181.0 




23.25+0.06 


22.20+0.01 




1.04+0.01 




19 


12:19:23.044 


+05:47:52.45 


1.59 


2.30 


0.98 


75.1 


25.39+0.04 


24.27+0.06 


23.06+0.01 


1.20+0.03 


1.26+0.02 




20 


12:19:23.052 


+05:47:53.10 


1.97 


1.89 


0.81 


124.6 


24.63+0.03 


23.67+0.06 


22.80+0.01 


1.05+0.02 


0.78+0.01 




21 


12:19:26.371 


+05:47:52.93 


1.07 


2.17 


0.96 


72.4 


26.03+0.06 


24.39+0.06 


23.16+0.01 


1.68+0.04 


1.29+0.02 




22 


12:19:26.542 


+05:47:53.69 


1.25 


2.44 


1.00 


147.1 


24.37+0.03 


23.53+0.06 


22.45+0.01 


0.93+0.01 


1.09+0.01 




23 


12:19:23.591 


+05:47:54.44 


1.05 


1.82 


0.95 


122.4 


25.12+0.04 


23.71+0.06 


22.51+0.01 


1.47+0.03 


1.24+0.01 




24 


12:19:25.284 


+05:47:54.66 


1.39 


2.38 


0.97 


85.1 


24.84+0.03 


24.17+0.06 


23.14+0.01 


0.77+0.02 


1.01+0.02 





(1) Globular Cluster ID 

(2) Equatorial Right Ascension (J2000) 
('^) Equatorial Declination (J2000) 

Elongation (major to minor axis ratio) 
(5) FWHM assuming a Gaussian core, in units of WFPG2 pixels 

SExTRACTOR stellarity index (1 = point like; = extended) 

(S/N) ratio within aperture 
(®) Aperture-corrected B band magnitude (if GC detected in B band) 

Aperture-corrected V band magnitude 
(-'^"^ Aperture-corrected / band magnitude 

(11) B-V color (if GC detected in B band) 

(12) v-l color 
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APPENDIX A: INCOMPLETENESS SIMULATION 

In order to evaluate the incompleteness of the GCs samples (i.e., the fraction of GCs not detected due to faintness 
or issues related to the detection process), we set up an artificial source test to calibrate the results from SEXTRACTOR. 
Simulated GCs were added to the NGC 4261 HST mosaics and their characteristics were measured with SEXTRACTOR using 
the same setup as for the real data. The simulation was repeated several times in order to improve the statistical results. In 
the following, we describe the simulation scheme: 

Generating GC radial profiles. We adopted the King model l|Kinelll963 . equation 14) to describe the GC brightness profile: 



where Jc is the scale brightness, rc is the core radius (roughly similar to the half light radius r^), r* is the tidal radius 
(which defines the radius up to which the gravitational field of the GC system can overcome the galactic gravitational field). 
The model can also be defined in terms of Vc and the concentration parameter c = logio{rt/rc), which expresses the core 
extent in terms of the GC size. In our analysis, we adopted r^ and c as the parameters characterizing the model. Image 
templates reproducing King profiles, as would be seen at the distance of NGC 4261, were created spanning the typical range 
of rc (1.5 < Tc < 10.0 pc) and c (1.0 < c < 2.0). We sampled from this 2-dimensional library of models when creating 
an artificial object. The c-space was sampled uniformly since catalogues of structural parameters for the GC popu lations 
in ne arby galaxies show that the probability density distri bution of c is almost fiat between the imposed limits (e.g. iHarrisI 
Il996l for the Milky Way: iBarmbv Holland, fc Huchra|[20o3 for M 31). We sampled the re-spac e following the ru distribution 
function measured by Jordan for the GC population of early type galaxies in the Virgo cluster l|jordan et al.li2005. . equations 
23 and 24). 

Generating instrument PSF. Artificial PSFs for different locations on each CCD of the HST cameras were produced using 
the Tiny Tim tool (jKristj i2995i) , assuming a black body spectral distribution peaked at 3500 K (close to the typical colour 
temperatures of stars composing a GC). In order to investigate the degree of variation of the PSFs across the WFPC2 field, we 
generated PSFs at the corners and the center of each CCD of WFPC2 and we measured their Encircled Energy as a function 
of radius. The aperture used for the photometry on the real data encompasses 95% of the PSF flux, with variations of only 
a few percent between different CCD locations or even different CCDs. Since these differences would not significantly affect 
the results of our simulation, we decided to adopt one single PSF model for the whole field, corresponding to the PSF at the 
central pixel of a WFC CCD. 

Defining source brightness. A set of sources of different brightness was generated sampling between 26.0 and 19.0 mag 
(these limits were chosen to cover the whole range of observed magnitudes in each band). The sampling was slightly finer for 
fainter magnitudes, where most deviation between input and output parameters is expected. 

Adding artificial sources on the NGG 4'i61 mosaics. For each iteration, a catalogue of uniformly distributed artificial 
source positions was generated. One random magnitude and one random King model from the set was assigned to each object 
position, without assuming any correlation between the size and the brightness of an artificial source. We convolved the King 
models with the instrument PSF and added the result on the NGC 4261 mosaic using the IRAF task mkobjects (IRAF package 
2.13-BETA2, 2006). The spatial density of artificial objects was limited to 1 object every 10 times the physical size of the 
simulated King model in order to avoid source confusion between the simulated objects themselves (each simulated object 
must be considered as an independent trial - many objects were simulated at the same time just to speed up the process). 
Simulated objects, though, could be cast onto real sources in the field, therefore accounting for real source confusion. The 
objects were uniformly distributed within the elliptical annulus constrained between the galactocentric radii Roc ~ 25" (GC 
detection limit imposed by the background) and Roc ~ R25 (~2'; galaxy major semi-diameter), as defined along the major 
axis. The axis ratio and P. A. parameters of the annulus have been obtained from a 2D model of NGC 4261 (see H4.1.2p . 

Artificial source photometry. Each simulated field was analysed by SEXTRACTOR in order to obtain the aperture fiux, 
position, axial ratio, FWHM and stellarity index of the artificial objects. The main aim of the simulation was to evaluate these 
parameters against the corresponding inputs as a function of their brightness. As a byproduct, we also got the incompleteness 
factors (i.e., the percentage of undetected sources) as a function of magnitude. 

Applying selection. The criteria used to define the "GC candidates" sample ( i]3.ip were applied to the SEXTRACTOR 
catalogues, in order to reproduce the same selection as for the real data. 

Statistics. For each iteration, we estimated the median of the measured values for each parameter, for the objects of 
a given input magnitude (the mean is not reliable since it is strongly driven by the outliers). At the end, the average of 
the medians of all the trials was selected as the expected value for the parameter. We estimated the dispersion around the 
expected value as the standard deviation of the medians. 

The results of the simulation are discussed in H3.2I and plotted in Figure [2] and O 
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Globular Cluster ID 
Equatorial Right Ascension (J2000) 
Equatorial Declination (J2000) 
'^^^ Elongation (major to minor axis ratio) 

FWHM assuming a Gaussian core, in units of WFPC2 pixels 
SExTRACTOR stellarity index (1 = point like; = extended) 
''^^ (S/N) ratio within aperture 

Aperture-corrected B band magnitude (if GC detected in B band) 
Aperture-corrected V band magnitude 

Aperture-corrected I band magnitude 

B-V color (if GC detected in B band) 

V-I color 



